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Abstract 

This communication presents a new approach to the re- 
construction of epicardial potentials from measured body 
surface potentials. Its salient characteristic is thai the 
temporal evolution of the epicardial potentials is ac- 
counted for by means of an explicit rttodel. The type of 
the model is chosen so as to yield efficient reconstruction 
using Kalman filtering. Classical Tikhonov regularization 
appears as a special case of the resulting technique, which 
can also be viewed as a generalization of the work pre- 
sented in [1]. Results obtained on real data and compari- 
son with those produced by existing methods confirm the 
interest of the approach. 


1 Introduction 

Our goal is to reconstruct epieardial potentials (EPs), 
whose temporal and spatial characteristics are of inter- 
est, from measured body surface potentials. Since the 
corresponding inverse problem is ill-posed, it must be reg- 
ularized through incorporation of prior information on the 
solution. Tikhonov regularization, which presents the ad- 
vantage of providing a solution in closed form, has been 
widely used [2|. However, in this approach, the time- 
continuity of the cardiac activation process is not ac- 
counted for as the data arc independently processed i hno- 
frame by time-frame. Oster &. Rudy [l] assessed the im- 
portance of temporal information and round it could im- 
prove the results significantly. However, practical use of 
temporal information remains an open question. In our 
iipproach. the time-dependence across time-frames is ac- 
counted for through a linear prediction model, which al- 
lows us to reconstruct the EPs using a Kalman filler 

2 Methods 

The relationship between body surface potentials and IiTs 
can be considered linear with a good approximation [3]. 
1 hcrefore, at time sample k, one can write 

*r(fc) = tf#jj(fc) I n(k), (1) 


where # r and <P E respectively denote the vectors of body 
surface potentials and EPs, and where H is the transfer 
matrix between the heart and the torso; n is a Gaus- 
sian white noise vector which represents modeling and 
measurement errors. In order to account for the time- 
con elation between KPs, the followii 
ear prediction model is introduced: 


#t-(fc I 1) = F# E (fc) + e(fc+l), 


(2) 


where F and e denote the time-invariant prediction ma- 
trix and the prediction error vector, respectively. (2) 
and (1) form a state-space representation of the phenom- 
ena, wt icli a/lows us to compute the least mean squares 
cstimavs of the KPs. with the following Kalman filter (ex- 
pressed here in compact form): 

*^(A- + 1 | k) •- Ff R (k | k - 1) 4 K k R k ' 

(#r(fc)-H#^(fc|*-l)) (3) 
K k $ FP klk ,H T (4) 
R k - IT + HP k}k X H T (5) 
P k .v>. -- FP k]k .- l F T - K k R k l K T k + R e .(6) 

In (3)-(6), R" and R e denote the covariance matrices of 
processes n and e. These matrices, along with the initial 
values of the estimate | 0) and of its covariance 

matrix P,| 0 , control the behavior of the Kalman filter. 

At this point, two major problems need to be solved: 
(0 determination of the stale-space model (2)-(l), i.e., 
H and F; (it) specification or the quantities which con- 
trol the Kalman filter. Regarding point (i), determina- 
tion of H has been widely studied. Mere, this matrix 
was computed using a 3-D finite element model or the 
torso according to the technique described in [3]. Deter- 
minntion of F is less classical, and we chose to estimate 
its value from recorded EPs. This estimation problem 
is also ill-posed, and therefore must be regularized. As 
the degree of accuracy of t tie prediction model necessary 
to achieve satisfactory reconstruction was difficult to as- 
sess beforehand, two types of regularizing constraints were 
used. I) F = nl, where n was estimated using a least- 
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squares method. 2) Determination of F using a regu- 
larized least-squares approach. The regularizing matrix 
was selected so as to make the predicted potentials at a 
given electrode depend mostly on its neighbors. The es- 
timator was implemented using a recursive least-squares 
algorithm. In order to solve point (ii), the prediction er- 
ror was analyzed so as to determine the corresponding 
covariance matrix R*. Setting the initial time just before 
cardiac activation allowed us to initialize the Kahnan fil 
ter with | 0) = 0 and P lj0 = al,a < 1. Finally, 

n was assumed uncorrelated spatially which gives R n a 
diagonal structure. Its entries were specified heuristically 
according to the observation signal-to-noise ratio (SNR). 

3 Results 

EPs were measured with a 03-lead mapping system using 
a sock electrode array in 8 patients with Wolff-Parkinson- 
White syndrome who had undergone arrhythmia surgery. 
Each recording was taken during normal sinus rhythm 
and lasted 1.024 s. In order to assess the generality of 
the prediction model, F was identified for each patient 
separately and used for predicting the EPs of the ot her 
patients, for both types matrices F. The results were 
compared using a relative squared error (USE) criterion. 

For general matrices F, estimation and cross-patient 
RSEs ranged from 0.01 to 0.25. and from 0.38 to 0.80, 
respectively. This indicates that F is specific to each pa- 
tient. For F = al, estimation and cross-patient RSEs 
were in the intervals [0.37,0.50j and [0.20,0.60] respec- 
tively, which shows that this prediction model is less ac- 
curate, but more widely usable than the previous one. 

Body surface potentials were simulated by multiplying 
EPs with transfer matrix H and adding Gaussian noise to 
the result. Two complete data sets were generated, with 
respective SNR values of infinity and 20 dD. Then re- 
construction of the EPs was carried out using Kalman fil- 
ter (3)-(6). Patient-specific prediction models were used, 
both for general matrices F and for F = al. R n was set 
as though the SNR value were 100 dB for the first data set, 
and 10 dB for the second one. For the sake of compari- 
son, two other reconstruction methods were implemented: 
classical Tikhonov regularization, and Twomey regular- 
ization using 4>fi(fc — 1) as the estimate of the solution [1] 
In these cases, the regularization parameter was deter- 
mined so as to minimize the reconstruction RSE. The 
results are collected in Table 1 . Best reconstructions were 
obtained using our approach with a general matrix F 

4 Discussion and conclusion 

The above results illustrate the interest of our approach, 
which produced low reconstruction errors. Clearly, the 
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critical point, is the selection of an appropriate predx 1 1< 
model. As F is patient-specific, the most satisfactory ;>;> 
proach would be to estimate its value directly from W.i, 
surface potentials. Techniques developed for blind ru- 
ination could be used for this purpose, but simplificatim. 
of the structure of F would certainly be required in ci- 
der to reduce the number of unknown parameters. Nuti 
that the same techniques could also be used for deter 
initiation of the diagonal entries of R" . Other potential 
improvements are the use of a Kalman smoother so as n> 
further reduce the reconstruction RSE, and utilization of 
fast or asymptotic Kalman ('liters in order to decrease tin- 
computational load. 
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